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In the last decade, thanks to the development of sophisticated numerical codes, major 
breakthroughs have been achieved in our understanding of the formation of asteroid families by 
catastrophic disruption of large parent bodies. In this review, we describe numerical simulations 
of asteroid collisions that reproduced the main properties of families, accounting for both the 
fragmentation of an asteroid at the time of impact and the subsequent gravitational interactions 
of the generated fragments. The simulations demonstrate that the catastrophic disruption of 
bodies larger than a few hundred meters in diameter leads to the formation of large aggregates 
due to gravitational reaccumulation of smaller fragments, which helps explain the presence 
of large members within asteroid families. Thus, for the first time, numerical simulations 
successfully reproduced the sizes and ejection velocities of members of representative families. 
Moreover, the simulations provide constraints on the family dynamical histories and on the 
possible internal structure of family members and their parent bodies. 


1. INTRODUCTION 

Observed asteroid families in the main asteroid belt are 
each composed of bodies that are thought to originate from 
the catastrophic disruption of larger parent bodies (e.g., 
Farinella et al., 1996). Cratering collisions can also lead 
to families, such as the one associated with asteroid Vesta, 
but we do not address this origin scenario here as few fami¬ 
lies have been linked to cratering events and their modeling 
requires a different approach (see the chapter by Jutzi et 
al. in this volume for more details). A few tens of aster¬ 
oid families have been identified, corresponding to groups 
of small bodies well-concentrated in proper-orbital-element 
space (see, e.g., Hirayama 1918, Arnold 1969, and the 
chapter by Nesvorny et al. in this volume) and sharing sim¬ 
ilar spectral properties (see, e.g.. Chapman et al. 1989 and 
the chapter by Masiero et al. in this volume). Large fami¬ 
lies contain up to several hundred identified members, while 
small and compact families have of the order of ten identi¬ 
fied members. Interestingly, the theory of the collisional 
origin of asteroid families rested for decades entirely on 
these similarities in dynamical and spectral properties and 
not on the detailed understanding of the collisional physics 


itself. Indeed, laboratory experiments on centimeter-scale 
targets, analytical scaling rules, or even complete numerical 
simulations of asteroid collisions were unable to reproduce 
the physical and dynamical properties of asteroid families 
(e.g., Ryan and Melosh, 1998). The extrapolation of labo¬ 
ratory experiments to asteroidal scales yielded bodies much 
too weak to account for both the size distribution and the 
dynamical properties of family members. In other words, 
there was no solution to match both the sizes and ejection 
velocities of family members simultaneously. To produce 
the large (assumed coherent) fragments seen in real fam¬ 
ilies required an impact energy leading to ejection speeds 
of individual fragments that were much too small for them 
to overcome their own gravitational attraction. The parent 
body would then be merely shattered but not dispersed and 
therefore no family would be created (Davis et al. 1979). 
Conversely, matching individual ejection velocities and de¬ 
riving the necessary fragment distribution resulted in a size 
distribution in which no big fragment was present, contrary 
to most real families (e.g., Davis et al., 1985, Chapman et 
al., 1989). Thus, big sizes implied no fragment dispersion 
(at the level of the dispersion of family members), and frag¬ 
ment dispersion implied no big fragments (at the level of 


1 


the sizes of family members). 

A big caveat in the extrapolation of laboratory results 
to large asteroid scales is that the role of gravity (of both 
the targets and their fragments) is not taken into account 
in a laboratory-scale disruption involving cm-size targets. 
Thus, the role of gravity in the catastrophic disruption of a 
large asteroid at the origin of a family remained to be as¬ 
sessed. Indeed, a possible scenario reconciling the sizes 
and ejection velocities of family members could be that the 
parent body (up to several hundred kilometers in size) is 
disrupted into small pieces by the propagation of cracks re¬ 
sulting from a hypervelocity impact but then the small frag¬ 
ments generated this way would typically escape from the 
parent and, due to their mutual gravitational attraction, reac¬ 
cumulate elsewhere in groups in order to build up the most 
massive family members. This idea had been suggested pre¬ 
viously by Chapman et al. (1982), and numerical simula¬ 
tions by Benz and Asphaug (1999) had already shown that 
at least the largest remnant of an asteroid disruption had to 
be a bound aggregate. However, the formation of a full fam¬ 
ily by reaccumulation of smaller fragments remained to be 
demonstrated. 

In the last decade, the formation of asteroid families was 
simulated explicitly for the first time accounting for the two 
phases of a large-scale disruption: the fragmentation phase 
and the gravitational reaccumulation phase. In these sim¬ 
ulations, the two phases are usually computed separately 
using a hybrid approach. This chapter reviews the major 
advances achieved since Asteroids III thanks to this new 
modeling work, and the implications for our understanding 
of asteroid family formation and properties. 

2. SIMULATING A FAMILY-FORMING EVENT 

Families are thought to form from the disruption of a 
large asteroid, called the parent body, as a result of the im¬ 
pact of a smaller projectile. Simulating such a process re¬ 
quires accounting for both the propagation of cracks in the 
parent body, leading to its conversion into separate frag¬ 
ments, and the possible gravitational interactions of these 
fragments. As explained above, the latter gravitational 
phase turns out to be crucial for reproducing asteroid fam¬ 
ily properties. In asteroid disruptions resulting from an hy¬ 
pervelocity impact, the fragmentation and the gravitational 
reaccumulation phases have very different associated dy¬ 
namical times. In the fragmentation phase, the time scale 
for the propagation of the shock wave is determined by the 
target’s diameter divided by the sound speed of the material 
(a few to tens of seconds for an asteroid 100 km in diam¬ 
eter). In the second phase, the time scale for gravitational 
reaccumulation is proportional to 1/> /Up, where G is the 
gravitational constant and p is the target bulk density, which 
corresponds to at least hours for p = 1 to 3 g cm -3 . There¬ 
fore it is possible to model the collisional event by separat¬ 
ing the two phases. A hybrid approach is generally adopted 
that consists of simulating first the fragmentation phase us¬ 
ing an appropriate fragmentation code (called a hydrocode; 


see the chapter by Jutzi et al. in this volume), and then the 
gravitational phase, during which the fragments produced 
by the fragmentation phase can interact under their mutual 
attraction, using a gravitational iV-body code. 

2.1 The Fragmentation Phase 

Several hydrocodes exist and are used in the planetary 
science community (see the chapters by Asphaug et al. and 
Jutzi et al. in this volume). In the first studies devoted to di¬ 
rect simulations of asteroid family formation ( Michel et al., 
2001, 2002, 2003, 2004), a three-dimensional “smoothed- 
particle hydrodynamics” (SPH) code was used. This code 
solves in a Lagrangian framework the usual conservation 
equations (mass, momentum, and energy) in which the 
stress tensor has a nondiagonal part. The first families mod¬ 
eled in this way (Eunomia, Koronis and Flora) were of S 
taxonomic type. S-type asteroids are expected to be mostly 
made of ordinary chondrite materials, for which basalt plau¬ 
sibly has similar properties and therefore the parent bod¬ 
ies of these S-type families were assumed to be non-porous 
basalt. The Tillotson equation of state for basalt was used 
(Tillotson , 1962), which is computationally expedient while 
sophisticated enough to allow its application over a wide 
range of physical conditions. Plasticity was introduced by 
modifying the stresses beyond the elastic limit with a von 
Mises yielding relation (Benz and Asphaug, 1994, 1995). 
A yielding relation accounting for the dependence of shear 
stress on pressure, such as the Mohr-Coulomb or Drucker 
Prager ones, is generally more appropriate for rock material 
(see the chapter by Jutzi et al. in this volume). It turns out to 
be important for cratering events for which part of the pro¬ 
cess is dominated by shearing or when the impacted body is 
composed of interacting boulders (a so-called “rubble pile” 
(Davis et al., 1979) or gravitational aggregate (Richard¬ 
son et al., 2002); see also Jutzi, 2015). However, it was 
found that in the case of the disruption of monolithic par¬ 
ent bodies, the details of the strength model (e.g., pressure- 
dependent vs. pressure-independent yield strength) do not 
play an important role (Jutzi 2015). For the lower ten¬ 
sile stresses associated with brittle failure, a fracture model 
was used, based on the nucleation of incipient flaws whose 
number density is given by a Weibull distribution (Weibull, 
1939, Jaeger and Cook, 1969). Durda et al. (2004) and 
Nesvorny et al. (2006) used a similar hydrocode to study the 
formation of satellites from asteroid disruptions and other 
family-forming events. Leinhardt and Stewart (2009) stud¬ 
ied large-scale disruptions and modeled the shock deforma¬ 
tion with an Eulerian shock-physics code, CTH (McGlaun 
et al., 1990), instead of the Lagrangian SPH code used in 
previous works. 

Recently, the SPH impact code used by Michel et al. 
(2001, 2002, 2003, 2004) was extended to include a model 
adapted for microporous materials (Jutzi et al., 2008, 2009, 
this volume). The formation of asteroid families formed 
from a microporous parent body, such as for dark- (carbona¬ 
ceous) type families, could thus also be investigated (Jutzi et 
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al., 2010, Michel et al., 2011). Another study looked at the 
case of rubble-pile parent bodies (containing macroporous 
voids; Benavidez et al., 2012). 

2.2 The Gravitational Phase 

Once the fragmentation phase is over and fractures cease 
to propagate (within the first simulated tens of seconds), the 
hydrodynamic simulations are stopped and intact fragments 
are identified. For impact energies typical of asteroid dis¬ 
ruptions and for targets with a diameter typically greater 
than 1 km, it was found that the bodies are totally shattered 
into fragments of mass equal to the mass resolution of the 
simulations. In the first simulations performed in this way, 
the numerical resolution was limited to a few 10 5 particles 
and corresponded to minimum boulder diameters of about 
1 to 4 km, for a parent body of a few hundred kilometers in 
diameter. Thanks to increased computer performance, it is 
now possible to perform simulations with up to several mil¬ 
lion particles. However, the gain in particle size resolution 
is not dramatic and simulations are still limited to minimum 
fragment diameters of a few hundred meters for target di¬ 
ameter of a few hundred kilometers. Reaching fragment di¬ 
ameters down to meters or less is beyond the capabilities of 
current and probably near-future technologies. Only when 
the target’s diameter is in the few hundred meters range can 
this minimum size be reached, but unfortunately, no aster¬ 
oid family can be identified involving a parent body of such 
a small size. 

Once identified in the simulation outcome, the fragments 
and their corresponding velocity distributions are then fed 
into a gravitational iV-body code, which computes the grav¬ 
itational evolution of the system from the handoff point to 
subsequent hours or days of simulated time. Because the 
number of fragments is up to a few 10 6 , and their gravita¬ 
tional interaction as well as their potential collisions need to 
be computed over long periods of time (up to several sim¬ 
ulated days), a very efficient /V-body code is required to 
compute the dynamics. The most appropriate code to tackle 
this problem, which is the only one used so far by various 
groups to simulate the outcome of the gravitational phase 
of a collision, is the code called pkdgrav (see Richardson et 
al., 2000 for the first application of this code to solar system 
problems). This parallel hierarchical tree code was devel¬ 
oped originally for cosmological studies. Essentially, the 
tree component of the code provides a convenient means 
of consolidating forces exerted by distant particles, reduc¬ 
ing the computational cost, with the tradeoff of introducing 
a slight force error (of order 1%) that does not affect the 
results appreciably since the dynamics are dominated by 
dissipative collisions. The parallel component divides the 
work evenly among available processors, adjusting the load 
at each timestep according to the amount of work done in 
the previous force calculation. The code uses a straightfor¬ 
ward second-order leapfrog scheme for the integration and 
computes gravity moments from tree cells to hexadecapole 
order. 


For the purpose of computing the gravitational phase of 
an asteroid disruption during which the generated fragments 
can interact and collide with each other, collisions are iden¬ 
tified at each step with a fast neighbor-search algorithm 
in pkdgrav. Once a collision occurs, because the relative 
speeds are small enough (of the order of meters per sec¬ 
ond), it is assumed that no further fragmentation takes place 
between components generated during the fragmentation 
phase. In fact, the simulations presented by Michel et al. 
(2001) assumed perfect sticking of colliding fragments and 
all colliding fragments were forced to merge into a single 
particle regardless of their relative velocities. This assump¬ 
tion is justified because the initial impact results in an over¬ 
all expanding cloud of fragments of relatively small individ¬ 
ual masses, down to the minimum fragment size imposed 
by the numerical resolution, and colliding fragments have 
typical relative speeds that are smaller than their individual 
escape speeds. Since the fragments in pkdgrav are repre¬ 
sented by spheres, when two spherical fragments reacumu- 
late, they are merged into a single spherical particle with the 
same momentum. The same assumption was used by Durda 
et al. (2004) and Nesvorny et al. (2006) in their studies of 
satellite formation and other family forming events. In a 
second and subsequent papers, Michel et al. (2002, 2003, 
2004) improved their treatment of fragment collisions by 
using a merging criterion based on relative speed and an¬ 
gular momentum. In this case, fragments are allowed to 
merge only if their relative speed is smaller than their mu¬ 
tual escape speed and the resulting spin of the merged frag¬ 
ment is smaller than the threshold value for rotational fis¬ 
sion (based on a simple prescription of a test particle re¬ 
maining on the equator of a sphere). Non-merging colli¬ 
sions are modeled as bounces between hard spheres whose 
post-collision velocities are determined by the amount of 
dissipation taking place during the collisions. The latter is 
computed in these simulations using coefficients of restitu¬ 
tion in the normal and tangential directions (see Richardson, 
1994, for details on this computation). Note that Durda et 
al. (2011, 2013) performed bouncing experiments between 
1-meter granite spheres as well as between cm-scale rocky 
spheres. These experiments gave a value for the normal co¬ 
efficient of restitution of « 0.8, although much lower val¬ 
ues are found with increasing roughness of contact surfaces. 
These results are particularly interesting because they are 
performed in an appropriate size regime (meter-sized bod¬ 
ies). However, bouncing in simulations occur at somewhat 
higher speeds (up to tens of m/s) than in those experiments, 
which may result in a decrease in the coefficient of restitu¬ 
tion due to the start of cracking and other energy dissipation 
processes. Moreover, although in our numerical modeling, 
perfect spheres are used, it is reasonable to account for ac¬ 
tual irregularities of fragments formed during the fragmen¬ 
tation phase to set the value of the coefficient of restitution. 
Since the values of these coefficients are poorly constrained, 
we usually arbitrarily set them equal to 0.5, meaning, for 
example, the rebound speed is set to half the impact speed. 

More recently, Richardson et al. (2009) enhanced the 


3 



collision handling in pkdgrav to preserve shape and spin in¬ 
formation of reaccumulated bodies in high-resolution sim¬ 
ulations of asteroid family formation. Instead of merging, 
fragments are able to stick on contact and optionally bounce 
or subsequently detach, depending on user-selectable pa¬ 
rameters that include for the first time several prescriptions 
for variable material strength/cohesion. As a result, the 
reaccumulated fragments can take a wide range of shapes 
and spin states, which can be compared with those ob¬ 
served. This comes with a cost in terms of computation 
time as several weeks to months are needed for one simu¬ 
lation using a few tens of current processors. This is the 
reason why this approach has so far only been used for par¬ 
ticular cases, such as modeling the formation of the asteroid 
Itokawa ( Michel and Richardson, 2013), and not systemat¬ 
ically for family formation investigations. 

And finally, we must note that the Soft-Sphere Discrete 
Element Method (SSDEM) has been introduced in pkd¬ 
grav (Schwartz et al., 2012, chapter by Murdoch et al. 
in this volume), which accounts more realistically for the 
contact forces between colliding/reaccumulating particles. 
This method should eventually replace the one developed 
by Richardson et al. (2009) to investigate the shape of reac¬ 
cumulated fragments as it avoids arbitrary particle sticking 
and rather let the reaccumulated particles evolve naturally 
towards the resulting equilibrium shape of the aggregate. 
However, solving for all contact forces between particles 
over the whole time scale of the gravitational phase, and 
covering a large enough parameter space (accounting for 
the uncertainty on the various friction coefficients) remains 
a computational challenge. Nevertheless, some collisional 
studies started using the SSDEM implementation in pkd¬ 
grav focusing on low-speed impact events. In effect, no 
fragmentation code was used for the impact phase, which is 
needed for impacts during which the sound speed of the ma¬ 
terial is reached. Thus, Ballouz et al. (2014a, b) used pkd¬ 
grav and SSDEM to simulate low speed impacts between 
rotating aggregates and to investigate the influence of the 
initial rotation of colliding bodies on the impact outcome 
and the sensitivity of some friction parameters. The number 
of particles was small enough (10 4 at most) that simulations 
could be performed within a reasonable computation time. 
In the case of a family formation, the outcome of hydrocode 
simulations consists in several hundred thousands to mil¬ 
lions of particles. Feeding them into the SSDEM version 
of pkdrgav requires another level of computer performance, 
although tests are under way. 

3. MODELING THE FAMILY PARENT BODIES 

Different possible internal structures have been consid¬ 
ered for the family parent bodies. Monolithic parent bodies 
composed of one material type with or without microporos¬ 
ity (meaning micropores in the solid rock; see the chap¬ 
ter by Jutzi et al. in this volume, for a definition of mi¬ 
croporosity) have been considered, as well as pre-shattered 
or rubble-pile parent bodies, with or without microporosity 


in the solid components. The assumed pre-shattered state 
could be seen as a natural consequence of the collisional 
evolution of main-belt asteroids. Indeed, several studies 
(see, e.g., Asphaug et al., 2002; Davis et al., 2002; Richard¬ 
son et al., 2002) have indicated that for any asteroid, colli¬ 
sions at high impact energies leading to a disruption occur 
with a smaller frequency than collisions at lower impact en¬ 
ergies leading to shattering effects only. Thus, in general, 
a typical asteroid gets battered over time until a major col¬ 
lision eventually disrupts it into smaller dispersed pieces. 
Consequently, since the formation of an asteroid family cor¬ 
responds to the ultimate disruptive event of a large object, 
the internal structure of this body before its disruption may 
be shattered by all the smaller collisional events that it has 
suffered over its lifetime in the belt, as suggested by Housen 
(2009) based on laboratory experiments and extrapolations 
using scaling laws. This would result in the presence of 
internal macroscopic damaged zones and/or voids. 

To model a pre-shattered target, Michel et al. (2003, 
2004) devised an algorithm that distributes a given num¬ 
ber of internal fragments of arbitrary shape and size within 
the volume of the parent body. The reason the internal frag¬ 
ments are given arbitrary shapes is that a network of frac¬ 
tures inside a parent body resulting from many uncorrelated 
small impacts is unlikely to yield spherical internal frag¬ 
ments whose sizes follow a well-defined power law. Then, 
void spaces are created by randomly removing a given num¬ 
ber of particles from the fractured set. Since there are vari¬ 
ous ways to define a pre-shattered internal structure, Michel 
et al. (2004) also built a model of a pre-shattered parent 
body in which large fragments are preferentially distributed 
near the center and smaller fragments are generated close to 
the surface. 

Another model, closer to the definition of a rubble pile 
was also built by Michel et al. (2003). In that case, spher¬ 
ical components whose sizes followed a specified power 
law distribution were distributed at random inside the parent 
body. Particles not belonging to one of these spherical com¬ 
ponents were removed to create void space and particles 
at the interface of two or more spherical components were 
assigned to fractures. Some simulations were performed 
using those two additional models and the collisional out¬ 
comes did not show any major qualitative difference com¬ 
pared to those obtained from the first pre-shattered model. 
Benavidez et al. (2012) constructed arbitrary rubble-pile 
targets by filling the interior of a 100-km-diameter spherical 
shell with an uneven distribution of solid basalt spheres hav¬ 
ing diameters between 8 km and 20 km. However, simula¬ 
tions performed so far using such rubble-pile parent bodies 
used a version of a SPH hydrocode with a strength model 
that did not allow the proper modeling of friction between 
the individual components of the rubble pile. As found by 
Jutzi (2015), the bodies in this case show a fluid-like behav¬ 
ior and are very (somewhat unrealistically) weak. There¬ 
fore in the following, we will only consider the results ob¬ 
tained for monolithic and pre-shattered bodies (as defined 
by Michel et al., 2003, 2004) either with or without micro- 
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porosity. 

4. REPRODUCING WELL-KNOWN FAMILIES 

For the first time, Michel et al. (2001) simulated entirely 
and successfully the formation of asteroid families from 
monolithic basalt-like parent bodies in two extreme regimes 
of impact energy leading to either a small or a large mass 
ratio of the largest remnant to the parent body Mi r /M p b. 
Two well-identified families were used for comparison with 
simulations: the Eunomia family, with a 284 km diameter 
parent body and Mi r /M p b « 0.67, was used to represent 
the barely disruptive regime, whereas the Koronis family, 
with a 119 km diameter parent body and Mi r /M p b « 0.04, 
represented the highly catastrophic one. Both families are 
of S taxonomic type, for which ordinary chondrites are the 
best meteorite analog, but basalt material is typically used 
as analog material in collisional studies. In these simula¬ 
tions, the collisional process was carried out to late times 
(typically several days), during which the gravitational in¬ 
teractions between the fragments could eventually lead to 
the formation of self-gravitating aggregates (as a represen¬ 
tation of rubble piles) far from the largest remnant. These 
first simulations assumed perfect sticking of reaccumulated 
fragments, regardless of relative speed and mass. This treat¬ 
ment was improved by Michel et al. (2002), allowing for 
the dissipation of kinetic energy in such collisions and ap¬ 
plying an energy-based merging criterion, as described pre¬ 
viously. This improved treatment did not change the con¬ 
clusion obtained with the more simplistic method because 
typical relative speeds between ejected fragments are most 
often below their mutual escape speed. Therefore, this new 
set of simulations confirmed the idea that the reaccumula¬ 
tion process is at the origin of large family members. Durda 
et al. (2007) and Benavidez et al. (2012) made a systematic 
study of collisional disruption of monolithic and rubble-pile 
basalt-like 100 km-diameter parent bodies, assuming per¬ 
fect sticking during reaccumulation, over a large range of 
impact conditions, and then re-scaled their results to com¬ 
pare with real families, showing again that the reaccumula¬ 
tion process is necessary to find any good solution. How¬ 
ever, a caution about extending results from the disruption 
of 100 km-diameter parent bodies to observed families that 
originated from parent bodies very different in size from 
100 km is in order and was acknowledged by Durda et al. 
(2007). 

4.1 The Size Distribution of Family Members 

The role of geometric constraints accounting for the fi¬ 
nite volume of the parent body in the production of fam¬ 
ily members was investigated by Tanga et al. (1999) and 
Campo Bagatin and Petit (2001). By filling the parent body 
with spherical ( Tanga et al. 1999) or irregular ( Campo 
Bagatin and Petit 2001) pieces, starting from the largest 
member, they were able to reproduce the size distribu¬ 
tion of some asteroid families to an encouraging level of 


agreement. However, these models do not incorporate any 
physics, nor do they take into account fragment reaccumu¬ 
lation, therefore they do not provide any explanation for 
how family members are formed nor any prediction for their 
internal properties. Moreover, ejection velocities are not ad¬ 
dressed by this approach. 

The first full numerical simulations of catastrophic dis¬ 
ruption and gravitational reaccumulation by Michel et al. 
(2001) assumed a monolithic structure of the parent body 
represented by a sphere with material properties of basalt 
(no internal porosity was considered). These simulations 
already reproduced qualitatively the main properties of real 
family member size distributions. However, when looking 
into more quantitative details, it was found that the cumula¬ 
tive size distribution of simulated fragments was character¬ 
ized systematically by a lack of intermediate-sized bodies 
and a very steep slope for the smaller ones (see an example 
in Fig. 1). Such characteristics are not always observed in 
the size distributions of real family members. In fact, for 
some families, the size distribution looks rather continuous. 
In their systematic study, Durda et al. (2007), using the 
same internal structures (monolithic, basalt) but consider¬ 
ing 100 km-diameter parent bodies only, found a larger va¬ 
riety of size distributions in terms of power-law slopes and 
discontinuities, depending on the considered impact condi¬ 
tions. 

Nevertheless, it is also known that the outcome of a col¬ 
lision is influenced by the initial internal structure of the 
parent body and that, depending on the initial structure, the 
fragment size distribution may be more or less continuous. 
In order to check this, Michel et al. (2003, 2004) mod¬ 
eled parent bodies with an internal structure composed of 
different zones of voids and fractures, as if they had first 
been shattered during their collisional history before under¬ 
going a major event leading to their disruption, as explained 
in Sec. 3. The simulations of the Eunomia and Koronis 
family formations were redone using pre-shattered parent 
bodies ( Michel et al., 2004), and the results were compared 
with those obtained using monolithic parent bodies. The 
best agreement was actually found with pre-shattered par¬ 
ent bodies. In particular, in the case of the Koronis family, 
an interesting result was obtained from these simulations, 
which may have important implications concerning the real 
family history. The size distribution obtained from the dis¬ 
ruption of a pre-shattered parent body contains four largest 
fragments of approximately the same size, as can be seen 
in Fig. 1. This peculiar characteristic is shared by the real 
family, and has been a source of debate as it was assumed 
that a single collisional event cannot produce such a prop¬ 
erty (see Michel et al., 2004, for a discussion). Moreover, 
the simulation using a monolithic parent body did not result 
in such a distribution. It was then demonstrated numeri¬ 
cally for the first time, by using a pre-shattered parent body, 
that these fragments can actually be produced by the origi¬ 
nal event, and therefore no subsequent mechanism needs to 
be invoked to form them, which would otherwise require a 
revision of the entire family history (Marzari et al., 1995). 
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Koronis Family: Cumulative Size Distribution 



Fig. 1.— Cumulative diameter distributions in log-log plots for the fragments of the simulated Koronis families. The plot on the left 
was obtained with a projectile colliding head-on, whereas an impact with an angle of incidence 8 equal to 45° gave rise to that on the 
right. Different symbols are used to distinguish between parent-body models. The plots also show the estimated sizes of the actual 
members down to the completeness limit (Tanga et al., 1999). Note that the simulations using a pre-shattered target reproduce the four 
nearly identical largest members. 


According to these results showing that even old families 
may well have originated from pre-shattered parent bodies, 
it was concluded that most large objects in the present-day 
asteroid belt may well be pre-shattered or self-gravitating 
aggregates/rubble piles. 

Constraints provided by the measured size distribution of 
family members can eliminate formation scenarios in nu¬ 
merical simulations. An interesting example is the Karin 
cluster, a small asteroid family identified by Nesvorny et 
al. (2002) that formed ~ 5.8 Myr ago in the outer main 
belt. The estimated size distribution for this family, when 
first identified, was fairly smooth and continuous over all 
sizes. In particular, there was no big gap between the size 
of the largest member of the family, (832) Karin, and that 
of the second-largest member, (4507) 1990 FV. Numerical 
simulations by Michel et al. (2003) indicated that the best 
match to the continuous size distribution was provided by 
the break-up of a pre-shattered or rubble-pile parent body. 
Simulations starting from a monolithic parent body, on the 
other hand, produced size distributions showing a large gap 
between the sizes of the largest and next-largest fragments. 
The finding that the parent body of the Karin cluster needed 
to be a rubble pile was actually consistent with its his¬ 
tory. Specifically, the parent asteroid of the Karin cluster 
is thought to have been produced by an early disruptive 
collision that created the much larger Koronis family some 
2-3 Gyr ago. According to the results of Koronis family- 
formation simulations, the parent asteroid of the Karin clus¬ 
ter should have been formed as a rubble pile from Koronis 
family debris. 

However, Nesvorny et al. (2006) later revised the def¬ 
inition of the Karin cluster. In particular, they found that 
the original second-largest identified member of the family, 
(4507) 1990 FV, is in fact a background asteroid with no 
relation whatsoever to the recent breakup at the origin of 


the Karin cluster. Once this body is removed from the clus¬ 
ter membership, a large gap opens between the size of the 
largest family member and smaller members, a distribution 
that is now best reproduced in simulation by starting with a 
monolithic parent body. This change in implication for the 
internal structure of the parent body shows the importance 
of having a reliable estimate of the actual size distribution 
of family members. However, in the case of the Karin clus¬ 
ter, this change is problematic because the parent body of 
the Karin cluster is expected to be a rubble pile, if it is an 
original fragment of the Koronis-forming event. A solu¬ 
tion proposed by Nesvorny et al. (2006) is that the Karin 
cluster parent body was really formed by reaccumulation of 
smaller fragments during the Koronis family formation, as 
found in numerical simulations, but then, it was somehow 
consolidated into a more coherent body by various possible 
processes (lithification of regolith filling the interior, etc.). 
Another possibility is that in simulations, we are missing 
cases in which large intact fragments are created, so that the 
Karin cluster parent body could really have been a mono¬ 
lithic body. In fact, this systematic absence of large intact 
fragments in asteroid disruption simulations is often men¬ 
tioned as a potential issue when discussing, for example, the 
internal structure of Eros, imaged by the NEAR/Shoemaker 
spacecraft (see Sec. 5.1). 

A model of fragmentation adapted for microporous bod¬ 
ies (see the chapter by Jutzi et al. in this volume, for a 
definition of microporosity) has recently been implemented 
into an SPH hydrocode and tested against experiments on 
pumice targets (Jutzi et al., 2008, 2009). It then became 
possible to simulate the formation of asteroid families from 
a microporous parent body. A microporous structure is 
assumed to be appropriate for parent bodies of dark taxo¬ 
nomic type or primitive asteroid families. In effect, several 
pieces of evidence point to the presence of a high degree of 
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porosity in asteroids belonging to the C-complex, such as 
the low bulk density (ss 1.3 g/cm 3 ) estimated for some of 
them, for instance the asteroid (253) Mathilde encountered 
by the NEAR Shoemaker spacecraft ( Yeomans et al., 1997), 
and as inferred from meteorite analysis (Britt et al., 2006). 
This model adapted for microporous bodies was used to re¬ 
produce the formation of the Veritas family, which is clas¬ 
sified as a dark type family whose members have spectral 
characteristics of low-albedo, primitive bodies, from C to 
D taxonomic types (Di Martino et al., 1997). This family 
is located in the outer main belt and is named after its ap¬ 
parent largest constituent, the asteroid (490) Veritas. The 
family age was estimated by two independent studies to be 
quite young, around 8 Myr (Nesvorny et al., 2003, Tsiga- 
nis et al., 2007). Therefore, current properties of the family 
may retain signatures of the catastrophic disruption event 
that formed it. Michel et al. (2011) investigated the for¬ 
mation of the Veritas family by numerical simulations of 
catastrophic disruption of a 140-km-diameter parent body, 
which was considered to be the size of the original fam¬ 
ily parent body, made of either porous or non-porous mate¬ 
rial. Pumice material properties were used for the porous 
body, while basalt material properties were used for the 
non-porous body. Not one of these simulations was able 
to produce satisfactorily the estimated size distribution of 
real family members. Previous studies devoted to either 
the dynamics or the spectral properties of the Veritas family 
treated (490) Veritas as a special object that may be discon¬ 
nected from the family. Simulations of the Veritas family 
formation were then performed representing the family with 
all members except Veritas itself. For that case, the parent 
body was smaller (112 km in diameter), and a remarkable 
match was found between the simulation outcome, using a 
porous parent body, and the real family. Both the size distri¬ 
bution and the velocity dispersion of the real reduced fam¬ 
ily were reproduced, while the disruption of a non-porous 
parent body did not reproduce the observed properties very 
well (see Fig. 2). This finding was consistent with the C 
spectral type of family members, which suggests that the 
parent body was porous and showed the importance of mod¬ 
eling the effect of porosity in the fragmentation process. It 
was then concluded that it is very likely that the asteroid 
(490) Veritas and probably several other small members do 
not belong to the family as originally defined, and that the 
definition of this family should be revised. This example 
shows how numerical modeling can better constrain the def¬ 
inition of (or the belonging to) an asteroid family, provided 
(i) that the fragmentation model used to simulate its forma¬ 
tion is consistent with the possible material properties of the 
parent body, and (ii) that the family is young enough that a 
direct comparison with the modeling is possible. 

4.2 The Ejection Velocity Distribution of Family Mem¬ 
bers 

In addition to fragment sizes, numerical simulations also 
provide the ejection velocities. In general, impact sim¬ 


ulations find that smaller fragments tend to have greater 
ejection speeds than larger ones. However, there is still 
a wide spread of values for fragments of a given mass, 
which makes it difficult to define a power-law relationship 
between fragment masses and speeds, such as the ones of¬ 
ten used in collisional evolution models (see, e.g. Davis et 
al., 2002). Figure 3 shows an example of this relation for 
a simulation reproducing the Eunomia family as a result of 
the disruption of a monolithic (basalt-like) parent body im¬ 
pacted at an impact angle of 45°. 

In the case of real asteroid families, however, the dis¬ 
persion of their members is characterized through their 
orbital proper elements, in particular their proper semi¬ 
major axis, eccentricity, and inclination. These elements 
have long been assumed to be essentially constants of mo¬ 
tion that remain practically unchanged over astronomically 
long timescales (e.g., Milani and Knezevic, 1994), although 
some perturbations have been found to be capable of mod¬ 
ifying them, as we will explain below. Thus, we do not 
have direct access to the ejection velocities of family mem¬ 
bers. Fortunately, ejection velocities can be converted into 
a dispersion in orbital elements through Gauss’ equations 
(. Zappald et al., 1996), provided that both the true anomaly 
and the argument of perihelion of the family parent body at 
the impact instant are known or assumed. For a given fam¬ 
ily, the estimated values of the barycenter semimajor axis, 
eccentricity, and inclination can be used with the Gauss for¬ 
mulae up to first order in eccentricity to compute for each 
member the distance of its orbital elements 5a, 5e and 51 
from the barycenter of the family: 

2 

- [(1 + e b cos fo)Vr + e b sin/oVn], 

na b 1- el 

V 1 - e b \ e b + 2 cos So + e b cos 2 fo . 

-- -—- - -Vr+sin/oVB > 

nab 1 + eb cos /o 

V 1 “ e b cos(uj + fo) 

-— Vw, 

nab 1 + eb cos fo 

(i) 

where Vt, Vr and Vw are the components of the ejection 
velocity in the along-track, radial, and out-of-plane direc¬ 
tions, respectively, n is the mean motion, fo is the true 
anomaly of the parent body at the instant of the break-up 
and io is its argument of perihelion. Since these last two an¬ 
gles are not known, their values must be assumed. Zappald 
et al. (1996) showed that the most sensitive angle is f 0 . 
Assuming different values of this angle changes the shape 
of the cluster containing the family members in orbital ele¬ 
ment space. In other words, it defines whether the break-up 
generates a family that is spread in semimajor axis, in ec¬ 
centricity, or in inclination. 

Thanks to this conversion, it is thus possible to assess 
the realism of a numerical simulation of a family formation 
by comparing the dispersion of family members and sim¬ 
ulated fragments in the same space. Unfortunately, other 
mechanisms exist that, depending on the age of the consid¬ 
ered family, can obscure the original dispersion of family 
members. In fact, once a family is created, its members are 
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Veritas Family: Cumulative Size Distribution 




Fig. 2.— Cumulative size distributions of fragments from the simulations of the disruption of a Veritas monolithic parent body, either 
non-porous (left) or porous (right). The impact angle is 0° (head-on) and the impact speed is 5 km/s. The distribution of the real family 
is also shown for comparison. In this case, the family consists of all members except Veritas itself, which reduces the size of the parent 
body to 112 km. The simulated time is about 11.6 days after the impact. 


Eunomia (Monolithic) 


Projectile: R = 38 km, 0=45° 



Fig. 3.— Fragment diameter D (normalized to that of the parent body D p i) vs. ejection speed in a log-log plot obtained from a 
monolithic Eunomia parent body simulation using a projectile impacting at an angle of incidence 9 = 45°. Only fragments with size 
above the resolution limit (i.e., those that underwent at least one reaccumulation event) are shown here. From Michel et al. (2004). 






subjected to various perturbations. In particular, high-order 
secular resonances, mean-motion resonances even involv¬ 
ing multiple planets (Morbidelli and Nesvorny , 1999), and 
the Yarkovsky thermal effect ( Farinella and Vokrouhliclcy , 
1999) have been shown to be capable of altering the proper 
elements. Therefore, while proper elements have been as¬ 
sumed conventionally to retain the memory of the disrup¬ 
tion outcome conditions, these later studies demonstrated 
that this is not necessarily true, even for the proper semi¬ 
major axis in the case when the asteroid is small enough 
that Yarkovsky drift is effective (see Bottke et al., 2002, 
Vokrouhlicky et al., this volume). Depending on how old 
the family is, the current proper elements of family mem¬ 
bers cannot be interpreted as reflecting their starting con¬ 
ditions; rather, they must be seen as a result of such secu¬ 
lar processes acting over time, whose effects are to cause a 
slow diffusion of family members in orbital-element space, 
starting from a smaller dispersion. The Koronis family is a 
good example showing these effects. The current distribu¬ 
tion of Koronis family members in proper-element space is 
quite spread and its shape suggests that it has been subjected 
to the Yarkovsky effect as well as to the effects of nearby 
secular resonances and mean-motion resonances. Bottke 
et al. (2001) computed the dynamical evolutions of 210 
simulated Koronis family members under the influence of 
the Yarkovsky effect and dynamical diffusion due to sev¬ 
eral resonances (namely, the 5:2 and 7:3 mean-motion res¬ 
onances with Jupiter, a secular resonance that involves the 
precession rate of the small body’s longitude of perihelion 
g and the fundamental frequencies of Jupiter g-, and Saturn 
go). The test family members were started with a disper¬ 
sion that is consistent with the ones obtained from impact 
simulations of Koronis family formation. They were inte¬ 
grated over 700 Myr, which is still shorter than the esti¬ 
mated age of the family (> 1 Gyr). However, this evo¬ 
lution showed that the current shape of the family cluster 
in proper-element space does not represent the original one 
from the collisional event but is well explained by its sub¬ 
sequent evolution. 

Fortunately, if a family is young enough, its dispersion 
can still be close to the original one resulting from the par¬ 
ent body break-up, and in that case, the comparison be¬ 
tween numerical simulations of family formation and actual 
family dispersion is straightforward. On the other hand, the 
degree of spreading observed now, together with the knowl¬ 
edge of the degree of dispersion resulting directly from the 
break-up by numerical simulations, can better constrain the 
age of the family, once the efficiency of the diffusion pro¬ 
cesses is well assessed. Nesvorny et al. (2002, 2003) iden¬ 
tified several asteroid families with formation ages smaller 
than 10 Myr. These families represent nearly the direct 
outcome of disruptive asteroid collisions, because the ob¬ 
served remnants of such recent break-ups have apparently 
suffered limited dynamical and collisional erosion ( Bottke 
et al., 2005). As already described in the previous section, 
the Karin cluster and the Veritas family belong to this group 
of such young families. 


Fig. 4 shows an example in which the dispersion of the 
actual Veritas family is compared with that of fragments 
from an impact simulation of Veritas family formation. The 
simulated dispersion matches the shape of the ellipses rep¬ 
resenting the real dispersion. This result is consistent with 
the expectation that the orbital extent of the family is not 
produced by post-diffusion processes, which gives some 
credibility to family formation simulations that reproduce 
both the size distribution and velocity dispersion of actual 
members. 

5. IMPLICATIONS 

5.1 Internal Structure of Asteroids 

According to numerical simulations of family formation, 
all fragments produced by the catastrophic disruption of a 
large asteroid (typically larger than 1 km in diameter, e.g. 
in the gravity regime) consist of self-gravitating aggregates, 
except the smallest ones. If this is correct, then most as¬ 
teroids of at least second generation should be rubble piles. 
We note that Campo Bagatin et al. (2001) ran a number of 
simulations of main belt collisional evolution to assess the 
size range where reaccumulated bodies should be expected 
to be abundant in the main asteroid belt. They found that 
this diameter range goes from about 10 to 100 km, but may 
extend to smaller or larger bodies, depending on the prevail¬ 
ing collisional response parameters, such as the strength of 
the material, the strength scaling law, the fraction of kinetic 
energy of the impact transferred to the fragments, and the 
reaccumulation model. 

The collisional lifetime of bodies larger than a few tens 
to hundreds of kilometers in diameter is longer than the age 
of the Solar System, suggesting that most bodies in that size 
range are likely to be primordial, while smaller bodies are 
probably collisional fragments (see, e.g., Bottke et al. 2005 
and the chapter by Bottke et al. in this volume). The ex¬ 
act size above which a body is more likely to be primor¬ 
dial is somewhat model-dependent. Binzel et al. (1989), 
from a study of lightcurves, suggested that this transition 
occurs at a diameter of « 125 km. However, as this is a 
statistical measure, some smaller asteroids may still be pri¬ 
mordial and some larger ones may have broken up in the 
past. In fact, due to the variability in possible interior start¬ 
ing compositions, and variations in the chaotic dynamics of 
accumulation, the size above which a body is more likely 
to be primordial is dependent on the specific formation sce¬ 
nario, as well as the compositions, masses and velocities 
involved. Thus, some asteroids smaller than 100 km may 
still be primordial, and some larger ones may have broken 
up catastrophically in the past. This is especially true if 
one goes back to the very earliest formation, in the first 
few million years, when considering hit-and-run collisions 
(see the chapter by Asphaug et al. in this volume). These 
may have completely disrupted some of the largest aster¬ 
oids, as projectiles, when they experienced grazing colli¬ 
sions into larger target embryos. This makes the internal 


9 



0.065 


Veritas Family, D > 8 km, w + /=i80°, /=30° 

Projectile: u impact = 5 km/s; Angle of Incidence = 0°: a = 3 AU; i=0° 




Fig. 4.— Distribution of fragments larger than 8 km from a simulation of disruption of a porous monolithic Veritas family parent body 
excluding Veritas itself as a result of a head-on impact of a projectile at 5 km/s with orbital semi-major axis a of 3 AU and inclination 
/ of 0°. The outcome is represented in the ae plane (left) and a-1 plane (right). The superimposed ellipse is an equivelocity curve for 
speed cutoff of 40 m/s, parent body true anomaly / = 30° and argument of perihelion co = 150°. This curve was defined by Tsiganis 
et al. (2007) as that closest to representing the dispersion of the real Veritas family in orbital element space. From Michel et al. (2011). 


structure of middle-sized asteroids one of the most impor¬ 
tant aspects of these bodies that can be determined by fu¬ 
ture space missions and observations, allowing us to test 
our interpretations based on theoretical collisional studies 
(see also Scheeres et al., this volume). 

During the past 4 billion years, catastrophic disruption 
has been the result of hypervelocity collisions. Bottke et al. 
(2005) estimate that about 20 asteroid families have formed 
from the breakup of parent bodies larger than 100 km di¬ 
ameter over the last 4 billion years. But several hundred 
asteroids currently exist in the 100 km size range, making 
it likely that most of these are original bodies. In this re¬ 
gard, asteroid 21 Lutetia, approximately 90 km diameter, is 
a scientifically important object, of which we have obtained 
a quick glimpse during Rosetta’s 2010 flyby (see Barucci et 
al., this volume). The relatively high measured mass (bulk 
density 3.1 g/cm 3 )led Weiss et al. (2012) to interpret Lute¬ 
tia as being a partly differentiated, impact-shattered, but 
largely intact parent body, covered in a predominately chon- 
dritic outer component. Other interpretations are of course 
possible. 

Assuming that the transition between primordial and 
second-generation bodies occurs at diameters about 100 
km, what about Eros, whose diameter is much below this 
threshold and therefore should be a fragment of a larger 
body? There is still a debate about the internal structure 
of this asteroid as the images of its surface can be explained 
by either a fractured (but solid/strength-dominated) struc¬ 
ture or a rubble pile (Asphaug 2009, chapter by Marchi et 
al. in this volume). However, if Eros is not a rubble pile, 
its formation as a fragment of a large asteroid would need 
a solution that is not yet found in numerical simulations of 
catastrophic disruptions. Another point of view could thus 
be that it is a monolithic body that has been shattered in 
place (e.g. Housen 2009; Buczkowski et al. 2008). In this 


case major impacts fracture it in place, introducing only 
modest increases to its porosity. This requires very low 
strain rate of expansion, e.g. a small elastic strain at frac¬ 
ture, which may be consistent with size-dependent relation¬ 
ships for brittle failure. So this is probably feasible to form 
a shattered monolith when a single monolithic body is im¬ 
pacted, but with relatively low energy compared to disrup¬ 
tion. But then one must ask where did the single monolithic 
body come from to begin with, and why has it not been 
subsequently fragmented and jumbled by slightly more en¬ 
ergetic collisions? The alternative is that the grooves have 
nothing to do with brittle failure, but are instead planes of 
granular failure. 

Thus, so far, the formation of a dispersed cloud of siz¬ 
able fragments (larger than a few hundreds of meters) sys¬ 
tematically requires that the parent body is first fragmented 
into small pieces, down to the resolution limit of simula¬ 
tions (a few hundreds of meters), and then that gravitational 
reaccumulation takes place to form larger final remnants. 
This is probably what happened for Itokawa, which appears 
to be a rubble pile ( Fujiwara et al., 2006). In fact, if the 
major blocks on Itokawa were intact monolithic, then they 
would give us a kind of lower size range of intact fragments 
produced from large impacts. Using the version of pkd- 
grav with enhanced collision handling to preserve shape 
and spin information of reaccumulated bodies (Richard¬ 
son et al., 2009), Michel and Richardson (2013) showed 
that the process of catastrophic disruption and gravitational 
reaccumulation can form fragments with similar shapes as 
Itokawa’s shape and can explain the presence of a large 
amount of boulders on the surface, as observed. Fig. 5 
shows the outcome of such a simulation. We note that in 
this kind of modeling, the shapes of the aggregates formed 
by the reaccumulation process are parameter dependent. In 
particular, if we change the assumed strength of the aggre- 
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gates or the bouncing coefficients (see Michel and Richard¬ 
son, 2013, for a definition of these parameters), the final 
shape may be different. For instance, a preliminary simu¬ 
lation using a lower strength leads to a final largest aggre¬ 
gate that is more spherical. Because of the lower assumed 
strength, reaccumulating aggregates break more easily as a 
result of tidal and rotational forces, and therefore the ob¬ 
ject produced by this reaccumulation has difficulty keep¬ 
ing its irregular shape and instead becomes more and more 
rounded. Further studies are required to determine whether 
this type of outcome has some interesting implications, and 
to assess the actual sensitivity of the final shapes of reaccu¬ 
mulated objects to the parameters. It may be that we can 
provide some rough constraints on some of the mechanical 
properties of asteroids whose shapes are known, based on 
the parameters required to form them using this model. An 
extensive set of simulations is planned for this purpose, that 
will require long runs with current computer power. 

The production of rubble goes up with size, because it 
gets harder and harder (with increasing gravity) to liberate 
mass to escape speed than to beat it into small fragments 
that eventually can reaccumulate. The implication is that if 
Itokawa is a rubble pile, then Eros should be even more so. 
Whether this is represents reality awaits direct seismic or 
internal-structure exploration (e.g., by radar tomography) 
of objects Eros-sized and larger (see also Scheeres et al., 
this volume). 

5.2 Compositions 

Originally families were only identified on the basis of 
dynamical considerations. Then, once spectral observations 
became available, it was found that the vast majority of 
those families identified by dynamics showed remarkable 
homogeneous spectral properties (see Masiero et al., this 
volume). So, homogeneity in terms of spectral properties 
seems to be the norm. However, when an object satisfies 
the distance criterion to be associated with a family, it is 
often considered as an interloper when its spectral proper¬ 
ties do not match. Therefore, the identification of family 
membership also relies on homogeneous spectral proper¬ 
ties and whether the homogeneity in spectral properties is 
a reality or an assumption is not clear yet. Such homo¬ 
geneity can only be explained if the family parent body was 
homogenous itself, so that when fragments reaccumulate 
during the reaccumulation phase, there’s no mixture of dif¬ 
ferent materials taking place. Alternatively, it may also be 
that the reaccumulation process does not mix different ma¬ 
terials that could be initially present in the parent body or 
mixes it so well that the outcome still looks homogeneous. 
Otherwise, if the parent body was heterogeneous in com¬ 
position and if some mixtures happened, then the resulting 
family would show a variety of spectral properties within 
its members. In fact, if re-accumulation is a random pro¬ 
cess, we expect the particles of a given large fragment to 
originate from uncorrelated regions within the parent body. 
In that case, if the parent body was heterogeneous in com¬ 


position, then the composition of reaccumulated fragments 
could be a mixture of various material. Conversely, if the 
initial velocity field imposed by the fragmentation process 
determines the re-accumulation phase, the particles belong¬ 
ing to the same fragment should originate from well-defined 
areas inside the parent body. In addition, the position and 
extent of these regions would provide indications about the 
mixing occurring as a result of the re-accumulation process. 

Michel et al. (2004) traced back, at least for some of 
the largest fragments, the original positions within the par¬ 
ent body of the particles that end up forming the aggre¬ 
gates during some family formations. As an example, they 
traced the particles belonging to the three largest fragments 
of their simulation of the Koronis family formation back 
to their original positions inside the parent body. Recall 
that the Koronis family was formed in a highly catastrophic 
event, as its largest member is estimated to contain only 
4% of the parent body’s mass. In such an event, the re¬ 
accumulation process lasts up to several days, much longer 
than for a barely disruptive event, and gives rise to many 
gravitational encounters. Therefore, this kind of event may 
well lose the memory of the initial velocity field. Never¬ 
theless, it was found that particles forming a large reac¬ 
cumulated fragment originate from well-clustered regions 
within the parent body. This indicates that re-accumulation 
is definitely not a random process. Interestingly, the posi¬ 
tion of the original region depends greatly on the internal 
properties of the parent body. The largest remnant of the 
pre-shattered model of the Koronis parent body involves 
particles that were initially located between the core and 
the region antipodal to the impact point. Conversely, in the 
monolithic parent body, those particles were initially much 
more clustered in the core region, with no particles originat¬ 
ing from the antipode. This difference also holds true for 
the next largest fragments. Nevertheless, these results in¬ 
dicate that the velocity field arising from the fragmentation 
phase has a major influence on the reaccumulation process. 
Particles that eventually belong to a given fragment origi¬ 
nate from the same region inside the parent body. However, 
this location (as well as its extent, which determines the 
degree of mixing of the fragments) depends also on the par¬ 
ent body’s internal properties in a complex way. Recently 
Michel et al. (2015) looked at the cases of parent bodies 
with internal structures that could represent large asteroids 
formed early in the Solar System history. Some results are 
shown in Jutzi et al. (this volume). They confirm that most 
particles in each reaccumulated fragments are sampled from 
the same original region within the parent body. However, 
they also found that the extent of the original region varies 
considerably depending on the internal structure of the par¬ 
ent body and seems to shrink with its solidity. 

As a conclusion, the spectral homogeneity within a fam¬ 
ily may represent the material homogeneity of the initial 
parent body. It may also be due to the way reaccumula¬ 
tion takes place. But in that case and if the parent body was 
heterogeneous, although each family member would still be 
homogeneous, we may expect different spectral properties 
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Fig. 5.— Snapshots of the reaccumulation process following the disruption of a 25 km diameter asteroid. From left to right: first instant 
at the end of the fragmentation phase when all fragments (white dots) are about 200 m in diameter; the ejection of those fragments a few 
seconds later; the first reaccumulations that occur because of the slow relative speed between some fragments, showing the formation 
of a few aggregates represented by different grey levels; the formation of the largest fragment of this disruption by reaccumulation of 
several aggregates into a single body; and the final largest fragment shown at two different instants: the boulders on its surface and its 
overall shape are reminiscent of Itokawa. Credit: Michel and Richardson, A & A, 554. LI, 2013, reproduced with permission ©ESO. 


from one member to the other, depending on which orig¬ 
inal region of the parent body it samples. The fact that 
most families do not show strong spectral variations be¬ 
tween family members, at least in the data of ground-based 
observations and except if this is imposed by the member¬ 
ship criterion, is consistent with the idea of homogeneity of 
the family parent body. 

6. CONCLUSIONS 

Our understanding of the collisional physics and our ac¬ 
count for gravity in large asteroid disruptions have allowed 
numerical simulations to reproduce successfully the forma¬ 
tion of asteroid families, in agreement with the idea that 
these families originate from the disruption of a large par¬ 
ent body. Simulation results indicate that asteroid family 
members are not just the product of the fragmentation of 
the parent body, leading to intact fragments, but rather the 
outcome of the subsequent gravitational phase of the event, 
which allows some of the intact fragments to reaccumu¬ 
late and form gravitational aggregates, or rubble piles. For 
all considered cases so far (family parent bodies of diam¬ 
eter typically larger than tens of kilometers), this outcome 
is systematic for fragments larger than a few hundred me¬ 
ters. Thus, according to simulations, since it is believed 
that most bodies smaller than 100 km in diameter originate 
from the disruption of a larger body, then they should be 
rubble piles or heavily shattered bodies, which is consistent 
with the low measured bulk densities for some of them and 
the finding by Campo Bagatin et al. (2001) based on main 
belt collisional evolution modeling. Therefore, exploring 
the origin of asteroid families unexpectedly led to a result 
that has great implications for the entire asteroid population 
and its history. 

Moreover, it was also found that the outcome of a dis¬ 
ruption is very sensitive to the original internal structure of 
the parent body, in particular the kind and amount of inter¬ 
nal porosity. Thus, the comparison between simulation out¬ 
comes for various kinds of parent-body structures (mono¬ 
lithic, pre-shattered, microporous, rubble pile) and real fam¬ 
ily properties can help to constrain the internal properties of 
the parent body of the considered families and in the family 
identification itself. For instance, it was found that the Ver¬ 


itas family is very well reproduced if the asteroid Veritas 
itself is excluded from its family, which was already rec¬ 
ognized as a possibility before disruption simulations were 
performed to model the formation of this family. 

Thanks to the improved sensitivity of observations, al¬ 
lowing us to reach smaller asteroid sizes, and to the tools 
developed to better define asteroid families, new asteroid 
families keep being identified, especially small and young 
ones. The latter, which have not been affected yet by dy¬ 
namical diffusion or post-collisional processes, are a good 
test for numerical simulations that must be able to repro¬ 
duce them as they are. Such an exercise, which has already 
been done successfully for some young families (e.g., Ver¬ 
itas, Karin) must keep being performed as a check for our 
numerical models. In particular, new fragmentation models 
are continuously being developed, accounting for various 
possible strength models and fragmentation modes. Once 
they are validated at small laboratory scales by comparison 
with impact experiments, they can be used at large scale 
(with an associated fV-body code) to reproduce young fam¬ 
ily properties, allowing us to increase the range of internal 
structures and fragmentation modes that can be considered 
for the parent body. This modeling work, calling for dif¬ 
ferent models, is crucial to better constrain the possible in¬ 
ternal structure of family parent bodies, to refine the defini¬ 
tion of a family, and to understand whether some families 
are formed from differentiated/heterogeneous parent bod¬ 
ies, despite their apparent (or assumed) homogeneities. As¬ 
teroid families are very important tracers of the entire aster¬ 
oid belt history and as we have already seen, their under¬ 
standing can have profound implications on our determina¬ 
tion of the physical properties of asteroids in general. 

More work is also required to check in which context 
large intact fragments can be produced in numerical simu¬ 
lations of large asteroid disruptions. Although there is no 
firm conclusion about Eros’ internal structure, the fact that 
it may be a shattered object only and not a rubble pile raises 
the issue of the formation of such large fragments in a col¬ 
lision. It may also be that the reaccumulation process is 
followed by internal processes that may consolidate boul¬ 
ders together. Such processes would eventually transform 
reaccumulated fragments into a coherent body. If this were 
the case, then reaccumulation would not necessarily imply 
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a rubble-pile structure. Simulations of the reaccumulation 
phase now include the possibility to account for the final 
possible shapes of reaccumulated fragments. This model¬ 
ing needs further improvement to increase its realism but it 
will be very difficult, if even possible, to achieve the level of 
complexity needed to model the internal processes that may 
consolidate boulders together. Asteroid internal processes 
are poorly understood and depend on too many parameters 
and unknowns. 

Space missions dedicated to direct measurement of inter¬ 
nal structures, and possibly to their response to an impact 
(e.g., by using a kinetic impactor), are thus crucial to im¬ 
prove our understanding of these important internal proper¬ 
ties of asteroids and to check our modeling of the collisional 
and internal processes. Moreover, sample return missions 
as well as visits/fly-bys of members of asteroid families 
would also provide detailed information on their physical 
properties and would allow us to check whether ground- 
based measurements wash out some important data regard¬ 
ing their composition and possible variations among mem¬ 
bers. Asteroid families and the collisional process, which is 
at the heart of family formation and evolutionary main belt 
history, rely on our efforts to combine complex models and 
space- ground-based measurements. 
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